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20.  ABSTRACT  (Continued). 

Thimble  Shoal  portion  of  the  navigation  channel  consists  predominantly  of  non- 
cohesive  material  while  the  sediment  along  the  Elizabeth  River  portion  of  the 
navigation  project  consists  primarily  of  clays  and  silts.  Sedimentation  in  a 
third  portion  of  the  overall  project,  referred  to  as  the  Atlantic  Ocean  Chan¬ 
nel,  was  evaluated  analytically  without  using  a  numerical  sediment  transport 
model . 

Based  on  sedimentation  results  from  the  Elizabeth  River  numerical  model, 
the  increase  in  shoaling  caused  by  channel  deepening  as  proposed  will  be 
23  percent.  The  distribution  of  shoaled  material  will  not  be  significantly 
altered,  other  than  a  slight  increase  in  skewness  toward  the  downstream  end. 

Based  on  sedimentation  results  from  the  Thimble  Shoal  numerical  model, 
the  increase  in  shoaling  caused  by  channel  deepening  as  proposed  will  be  about 
20  percent.  The  distribution  of  shoaled  material  will  be  slightly  altered  in 
that  both  the  upper  and  lower  channel  shoaling  peaks  which  presently  exist 
will  tend  to  migrate  even  more  toward  the  ends  of  the  dredged  channel. 

Based  on  the  analytic  analysis,  the  estimate  of  shoaling  for  the  new 
Atlantic  Ocean  Channel  is  about  200,000  cu  yd  annually. 
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PREFACE 


In  December  1981,  <.he  US  Army  Engineer  Waterways  Experiment  Station 
(WES)  was  requested  by  the  US  Army  Engineer  District,  Norfolk,  to  conduct  an 
investigation  of  possible  sedimentation  changes  in  Norfolk  Harbor  and  Channels 
caused  by  the  proposed  deepening. 

The  study  was  conducted  by  personnel  of  the  Hydraulics  Laboratory,  WES, 
under  the  general  direction  of  Mess  s.  H.  B.  Simmons,  retired  former  Chief  of 
the  Hydraulics  Laboratory,  F.  A.  Herrmann,  Jr.,  Chief  of  the  Hydraulics  Labora¬ 
tory,  R.  A.  Sager,  Chief  of  the  Estuaries  Division,  E.  C.  McNair,  Chief  of  the 
Sedimentation  Branch,  and  R.  A.  Boland,  Chief  of  the  Hydrodynamics  Branch.  The 
project  was  conducted  and  this  report  prepared  by  Messrs.  R.  C.  Berger,  Jr., 

S.  B.  Heltzel,  R.  F.  Athow,  Jr.,  D.  R.  Richards,  and  M.  J.  Trawle.  Other  WES 
personnel  participating  in  the  study  were  Messrs.  J.  A.  Boyd,  D.  M.  Marzette, 
and  Ms.  V.  P.  Pankow. 

We  gratefully  acknowledge  the  valuable  contributions  of  Dr.  J,  C. 

Ludwick  of  Old  Dominion  University  and  Mr.  J.  R.  Melchor,  US  Army  Engineer 
District,  Norfolk. 

Commanders  and  Directors  of  WES  during  the  conduct  of  this  study  and 
preparation  and  publication  of  this  report  were  COL  Tilford  C.  Creel,  CE,  and 
COL  Robert  C.  Lee,  CE.  Technical  Director  was  Mr.  F.  R.  Brown. 
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CONVERSION  FACTORS,  US  CUSTOMARY  TO  METRIC  (SI) 
UNITS  OF  MEASUREMENT 


US  customary  units  of  measurement  used  in  this  report  can  be  converted  to 
metric  (SI)  units  as  follows: 


cubic  feet  per  second 
feet 


feet  per  second 
inches 

knots  (international) 
miles  (US  nautical) 
miles  (US  statute) 


0.4047 

0.02831685 

0.3048 

0.3048 

25.4 

0.514444 

1.852 

1.609344 


_ To  Obtain _ 

hectares 

cubic  metres  per  second 
metres 

metres  per  second 
millimetres 
metres  per  second 
kilometres 
kilometres 
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NORFOLK  HARBOR  AND  CHANNELS  DEEPENING  STUDY 


SEDIMENTATION  INVESTIGATION 
Chesapeake  Bay  Hydraulic  Model  Investigation 


PART  I:  INTRODUCTION 


Chesapeake  Bay 


1.  Chesapeake  Bay  with  its  tributary  estuaries  forms  the  largest  estua¬ 
rine  system  in  North  America.  The  190-mile*-long  estuary  varies  in  width  from 
4  to  30  miles  with  an  average  depth  of  28  ft  (Figure  1).  The  mean  annual  dis¬ 
charge  of  its  126  freshwater  tributaries  is  approximately  70,000  cfs,  almost 
90  percent  of  which  is  contributed  by  the  Susquehanna,  Potomac,  Rappahannock, 
York,  and  James  River  basins.  The  Atlantic  Ocean  provides  salt  water  to  the 
bay,  producing  large  salinity  variations  within  its  boundaries.  The  eastern 
shore  is  generally  saltier  than  the  western  shore,  attributed  in  part  to  the 
dominance  of  freshwater  flow  from  the  western  shore  tributaries  and  to  the 
counterclockwise  tendency  of  flow  resulting  from  Coriolis  force  (Richards  and 
Morton  1983). 

2.  Chesapeake  Bay  is  classified  geologically  as  a  drowned  river  valley 
estuary.  The  Holocene  sea-level  rise  inundated  the  Susquehanna  River  Valley 
to  form  the  bay.  Sedimentation  from  the  tributaries  as  well  as  erosion  of  the 
banks  has  contributed  to  maintaining  the  bay's  broad,  shallow  character.  The 
bay  is  classified  as  a  partially  mixed  estuary,  although  various  stages  of 
freshwater  discharge  and  tidal  and  wind  mixing  cause  portions  to  alternate  be¬ 
tween  well  mixed  and  highly  stratified.  Tides  are  semidiurnal  with  mean  ranges 
from  1  to  2  ft.  The  length  of  Chesapeake  Bay  is  such  that  a  complete  tidal 
wave  is  contained  within  its  limits  at  all  times.  Wind-generated  waves  are 
generally  less  than  3  ft  in  height,  but  larger  waves  can  occur  during  high 
wind  conditions.  Average  maximum  velocities  for  tide  and  wind-driven  currents 
range  from  0.5  to  3  fps  (Richards  and  Morton  1983). 

3.  Chesapeake  Bay  and  adjoining  tributaries  are  gradually  filling  with 


*  A  table  of  factors  for  converting  US  customary  units  of  measurement  to 
metric  (SI)  units  is  presented  on  page  3. 
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sediment.  The  sources  and  losses  of  suspended  material  within  this  estuarine 
system  are  complex  and  varied.  The  sources  of  this  sediment  include  the 
Susquehanna  River,  shore  erosion,  the  ocean,  and  biological  production.  In  the 
upper  reaches  of  the  bay,  upland  discharge  is  the  principal  source  of  sediment. 
However,  shore  erosion  is  the  major  source  of  sediment  in  the  middle  and  lower 
reaches  of  the  estuary.  These  processes  may  operate  according  to  a  seasonal 
cycle,  so  that  at  any  time  during  the  year,  one  of  these  processes  may  domi¬ 
nate  the  others.  For  instance,  the  Susquehanna  River  discharges  70  percent 
of  its  total  annual  sediment  load  during  the  spring  freshet.  In  some  cases, 
during  heavy  spring  floods,  the  river  will  discharge  the  bulk  of  its  sediment 
load  in  a  few  weeks.  Analysis  of  sediment  near  the  bay  mouth  indicates  mainly 
sandy  sediment.  Sediments  from  tributary  rivers  (i.e.,  James  River)  are  pre¬ 
dominant  silt  and  clay  with  progressively  coarser  sediment  toward  the  bay. 

The  tributaries  contribute  little  if  any  sand.  Most  of  the  coarse  sediment  is 
either  eroded  from  the  margins  of  the  lower  and  middle  portions  of  the  bay  or 
possibly,  according  to  some  opinions,  transferred  from  the  ocean  to  the  bay. 
Overall  hay  deposition  averages  0.8  mm  per  year. 

Norfolk  Harbor 

4.  Navigational  uses  of  Chesapeake  Ray  in  the  Norfolk  area  are  of  great 
importance  to  the  Nation  and  the  local  communities.  Due  to  its  naturally  pro- 
ti'ct.ed  harbors,  the  Norfolk  area  has  historically  been  the  home  port  of  naval 
activities  since  colonial  times.  Commercially,  Norfolk  has  played  a  major 
roio  in  east  coast  bulk  shipping  for  many  years.  Its  closeness  to  the  Ap- 
pn in chi  an  coal  fields  and  connecting  rail  lines  has  helped  it  become  the 
largest  coal  exporting  port  in  the  United  States.  However,  with  the  current 
trends  toward  deeper  draft  bulk  cargo  vessels  and  an  ever-increasing  demand 
for  United  States  coal ,  Norfolk  may  lose  some  of  this  competitive  advantage. 
Currently,  several  vessels  calling  on  Norfolk  must  carry  partial  loads  to  navi¬ 
gate  through  the  existing  channels.  Since  the  majority  of  the  cargo  passing 
through  Norfolk  is  high  in  volume  and  low  in  price,  the  efficient  use  of  ship¬ 
ping  is  crucial  to  bring  profits.  Unless  the  harbor  is  deepened,  future  deep- 
draft  vessels  may  he  forced  to  use  other  ports  (Richards  and  Morton  1983). 

Proposed  Channel  Improvements 

f>.  The  proposed  improvements  to  channels  and  anchorages  approaching 


Norfolk  Harbor  are  shown  in  Figure  2  and  described  as  follows  (USA!'.!,',  Norfolk, 
1980)  : 

a.  Increasing  the  depth  of  Thimble  Shoal  Channel  from  45  to  :  5  ft 

below  mean  low  water  over  its  existing  1,000- ft  width. 

h.  Increasing  the  depth  of  Norfolk  Harbor  Channel  from  45  to  3  5  ft 

below  mean  low  water  over  its  existing  800-  to  1  , 500  —ft  width  to 

the  coal  terminal  at  Lamberts  Point. 

c.  Increasing  the  depth  of  the  Channel  to  Newport  News  from  45  to 
55  ft  below  mean  low  water  over  its  existing  800-ft  width  to  the 
coal  terminal  at  Newport  News. 

d.  Dredging  a  new  channel,  referred  to  as  the  Atlantic  Ocean  Chan¬ 
nel,  off  Virginia  Beach  to  a  depth  of  57  ft  below  mean  low  water 
and  a  width  of  1,000  ft  over  a  length  of  10.6  miles. 

e.  Constructing  three  fixed -mooring  anchorage  facilities,  each 
capable  of  accommodating  two  large  vessels  simultaneously. 

_f.  Increasing  the  depth  of  the  Flizabeth  River  and  the  Southern 
Branch  of  the  Elizabeth  River  between  Lamberts  Point  (river 
mile  9)  and  the  Norfolk  and  Western  Railway  Bridge  (river 
mile  15)  from  40  to  45  ft  below  mean  low  water  over  Its  existing 
375-  to  750-ft  width. 

g.  Increasing  the  depth  of  the  Southern  Branch  of  the  Flizabeth 
River  between  the  Norfolk  and  Western  Railway  Bridge  (river 
mile  15)  and  the  US  Routes  460  and  13  highways  crossing  (river 
mile  17.5)  from  35  to  40  ft  below  mean  low  water  over  its  exist¬ 
ing  250-  to  500-ft  width,  and  providing  a  new  800-ft  turning 
basin  at  the  terminus  of  the  channel  improvement. 

6.  The  depths  listed  in  paragraph  5  are  project  depths  and  do  not  pro¬ 
vide  fur  advance  maintenance  or  dredging  tolerance.  The  actual  depths  for  the 
:  reposed  new  channels  with  the  advance  maintenance  and  dredging  tolerance 
should  be  3  ft  deeper.  Previous  deepening  projects  in  the  vicinity  also  had 
:r  ■  >v  i  s  i  eu for  similar  amounts  of  advance  maintenance  and  dredging  tolerance. 
The  Mimcrlr  tl  model  study  of  the  proposed  channel  deepening  used  the  existing 
urn  jci't  channel  depths  plus  an  additional  3  ft  to  allow  for  advance  maintenance 
and  dredging  tolerance  as  the  base  condition;  the  plan  condition  used  the 
deepened  pf"  ;cri  depths  plus  3  ft  for  advance  maintenance  and  dredging  toler¬ 
ance  (Richards  and  Morton  1983). 


Pn  rpose 

7.  The  purpose  of  this  study  wo.  to  investigate  the  impact  of  *  In  pj 
posed  channel  deepening  on  *  he  sol  'up  fa  Ion  rharari eri «l  i<  s  within  the 


Figure  2.  Project  map  showing  proposed  Norfolk  dredged  channels 


ication  and  node  num 


mesh  were  modeled  in  the  following  manner.  The  upstream  James  River  limit  was 
a  velocity  or  discharge  boundary,  as  was  the  Elizabeth  River  Western  Branch. 

The  southern  portion  of  Chesapeake  Bay  was  an  elevation  or  tidal  boundary.  The 
final  water  boundaries  were  the  locks  on  the  South  Branches  of  the  Elizabeth 
River,  represented  as  a  zero-velocity  boundary. 

31.  Report  1  of  this  series  discussed  the  results  of  the  physical  model 
testing  program.  Four  test  conditions  were  developed  as  representative  of 
conditions  commonly  occurring  in  the  prototype.  Velocity  and  tidal  elevations 
were  recorded  over  a  25-hr  period  for  each  test  condition  for  use  in  the  nu¬ 
merical  simulations.  The  test  conditions  are  summarized  below: 

Discharge,  cfs  Tidal  Range,  ft 

200,000  4.8 

200,000  3.0 

70,000  4.8 

70,000  3.0 

The  discharge  shown  above  was  the*  total  freshwater  discharge  into  the  Chesa¬ 
peake  Bay  fron  all  tributaries.  The  tidal  range  was  measured  at  the  Atlantic 
Ocean  physical  model  control  gage.  Figures  5  and  6  show  the  locations  of  the 
physical  model  velocity  and  water — surface  elevation  stations  and  their  respec¬ 
tive'  nodal  mesh  designations. 

32.  Initial  calibration  of  the  numerical  model  consisted  of  comparison 
of  Rase  4  'existing  channel  conditions  with  the  boundary  conditions  of  condi- 
t  Ion  *  )  numerical  water — surface  elevations  and  velocities  with  the  data  mea¬ 
sured  in  the  physical  model  at  the  same1  locations.  If  large  differences  were 
‘  usd  S- 1  we*")  th''  numerical  and  physical  model  results,  t  he  numerical  model 

p  trumet  ers  of  either  friction  or  eddy  diffusivity  were  adjusted  and  the  ihi- 
:■■»  ri-al  model  rerun.  After  satisfactory  agreement  had  been  achieved  at  all 
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PART  III:  MODELING  PROCEDURES 


Elizabeth  River 


Hydrodynamic  simulation 

27.  The  first  step  in  the  numerical  hydrodynamic  simulation  of  the  pro¬ 
totype  was  the  development  of  a  representative  grid  or  mesh.  Secondly,  the 
physical  model  data  were  examined  for  completeness  and  consistency,  and  then 
the  houndary  condition  files  were  generated  for  use  in  the  time-stepping  simu¬ 
lation.  The  final  steps  included  the  calibration  of  the  numerical  hydrody¬ 
namics  by  means  of  comparison  and  adjustment  for  best  fit  to  the  physical 
model  stations  in  the  interior  portions  of  the  numerical  mesh,  and  the  simu¬ 
lation  of  the  deepened  channel  configuration  for  all  boundary  conditions. 

28.  Mesh  development  incorporated  the  physical  boundaries  of  the  estu¬ 
ary  with  a  representation  of  the  bathymetry.  A  mesh  overlay  was  drawn  on 

the  most  recent  NOS  charts  available.  The  following  table  shows  the  charts 
used  and  their  respective  scales. 


NOS  Chart  No. 

Scale 

Date 

1  2253 

1:20,000 

Mar 

79 

12248 

1:40,000 

Nov 

74 

1  2222 

1:40,000 

Apr 

81 

The  mesh  intersections  and  midside  connections  between  intersections  were 
then  numbered  as  nodes,  and  the  mesh  inclosures  were  numbered  as  elements. 

The  Elizabeth  River  mesh  was  constructed  of  1,496  computational  nodes  and 
407  elements.  Figure  4  shows  the  resulting  mesh. 

29.  The  physical  model  data  were  supplied  to  the  authors  in  the  form 
of  magnetic  tape  and  a  tabular  printout.  The  tape  was  read  into  the  main¬ 
frame  computer  system  used  during  the  study  (Boeing  Computer  Services)  and 
the  format  was  changed  as  required.  A  preprocessing  code  (CODE  24)  was 
then  used  to  average  the  point-depth  values  to  yield  depth-integrated  values 
needed  for  the  hydrodynamic  code.  The  resulting  files  were  labeled  "boundary 
files"  and  were  used  as  boundary  updates  for  the  time-stepping  hydrodynamic 
simulations. 

10.  The  upstream  and  downstream  water  boundaries  of  the  Elizabeth  River 

1  7 


=  erosion  rate  constant 


.lit iiil  shear  stress  tor  particle  erosion 

C 

SITiJli  is  an  outgrowth  ot  the  mode  1  SEDIMENT  LI  (Ariathurai , 
Art  hur ,  .n.l  Krone  1  v>  7  7 )  developed  under  the  direction  of  R.  IS.  Krone  at 
.tiivers  it  ot  tali  tornia,  Davis. 


Mae- 

the 


22.  The  value  ol  Ceq  can  be  determined  i  ruin  an,  a;  s«.  vu ra *  t  r  r  : 

relations.  The  sand  version  of  the  sediment  model,  S'll.'iHl,  e  •  :  t  h.  r 
White  formula  (1973),  which  performed  sat  is  fact  or ;  lv  in  test  s  h  v  Wh.  ....!  ■>{  her 
(White,  Mi  1 1  i ,  and  Crabbe  1975  ,  Swart  197b).  The  t  ransport  potent  •  i  ’  • 
lated  to  sediment  and  flow  parameters  by  the  fol  lowing  expi  i  <to.  .  ' 
Ackers-White  formu  las . 

23.  Transport  potential  is  converted  to  t  ran*. port  arai  i  t  y  ’  ;  ai  :  ,  - 

ing  the  availability  constraint 

GS.  =  CP.  •  PI .  15) 

i  l  i 

where 

GS  =  transport  capacity 
i  =  number  of  grain-size  class 
GP  =  transport  potential 

PI  =  1/100  percent  at  bed  surface  covered  bv  grain  size,  class  i  ( i he 
percent  parameter  varies  from  ! 00  to  0  and  is  based  on  weight  c 
sediment  as  it  is  compatible  wirti  multiple  nr.  in-si /.e  clans 
concepts) 

24.  Clay  transport.  Deposition  rates  for  clay  beds  were  calculated 
with  the  equations  of  Krone  (1962). 


where 

Vs  =  fall  velocity  of  a  sand  grain 
D  =  f low  depth 

C  =  concentration  of  sediment 

T,  =  bed  shear  stress 
b 

’■  =  critical  shear  stress  for  deposition 

C  =  critical  concentration  =  300  mn/> 
c  '  /a 

Vk  =  Vs/C‘/J 
c 

25.  Krosion  rates  were  comput'd  by  a  s i mp 1 ; 1 i c at i on  of  Par t hen i a.h 
(1962)  for  part i c 1 e-by  •  part i e 1 e  erosion. 


Ih'  sou  ice  'or  i  s  i  ompai  <ul  by 


i  1  ov.'  vi'loriLv  in 


z-d  i  rect  ion  ,  m/sec 
direction  perpend  i  cu  1  ar  to  x  ,  m 

l 

c :  :  i->  •  t  i  ve  diffusion  coefficient  in  x-d i rect i on ,  in4" /hoc 
>'i  ii'i't  i'.'o  dit fusion  coefficient  in  z-d  i  root  ion 
a  coelt  ieient  for  tho  source  terra,  1/sec 

Lac  f.|iii  1  iliriura  concent  rat  ion  portion  of  the  source  term,  kg/m^/sec 
bed  stiear  stress.  I  lie  bed  shear  stress,  r  ,  takes  the  form: 

rh  =  pu*  (2) 


-  water  densiLy 
u ...  -  shear  velocity 

20.  The  Manning  shear  stress  equation.  The  Manning  form  of  the  shear 
sti'esi  equation  was  used  in  this  study 


=  un_ 

"  (IMF.  Dl/6 


(3) 


v,':. e  re 


u  flow  velocity 


a  =  Mann  i  in;' s  rouvimess  value 

0MK  =  coefficient  ot  1  for  metric  units  and  1.486  for  English  units 
acceleration  due  to  gravity 
J  ~  flow  depth 


Sami  transport.  the  supply  of  sediment  to  and  from  the  bed 
bed  outer  ini  (sand)  was  controlled  by  the  transport  potential 
i be  avuilabil i tv  ol  material  in  the  bed .  The  bed  source  term 


for  non- 
of  the 
is 


S 


Ccq  -  C 
t 

c 


(4) 


S  =  Si  ur.'e  I  ■ '  I  ". 

f  >  • ' ,  ~  '  lilihrinr,  concentration 
(  -  i  on.  ent  rat  i  si  of  sediment 

t  :  '  f  i r.K  t ei  : et  j  i  rime  tor  effecting  the  transition 


Shoal  study.  1'he  cohesive  or  clay  version  of  S'iTDH  was  used  in  the  Elizabeth 
River  study. 

Hie  jiydrodyn_amic  _mod_e_l_,  JtMAydY 

15.  The  hydrodynamic  model,  RMA-2V,  solves  the  depth- integrated  equa¬ 
tions  o!  conservation  of  mass  and  momentum  in  two  horizontal  directions.  Tlio 
present  model  is  an  improvement  of  an  earlier  version,  RMA-2  (Norton  and  King 

1177).  The  model  is  formulated  in  terms  of  velocities  and  turbulent  exchange 

coe  f  f i c ients  . 

Lb.  The  1  inite  element  method  using  C.alerkin  weighted  residuals  is  used 
to  solve  the  conservation  of  mass  and  momentum  equations.  Individual  elements 
may  be  either  quadrilaterals  or  triangles  and  may  have  curved  (parabolic) 

sides.  file  shape  functions  are  quadratic  for  flow  and  linear  for  depth.  In¬ 

tegration  in  space  is  Caussian.  Derivatives  in  time  are  replaced  by  nonlinear 
finite  difference  approximations. 

17.  The  finite  element  solution  is  fully  implicit  and  the  set  of  simul¬ 
taneous  equations  is  solved  by  Newton-Raphson  iteration.  The  solution  is 
achieved  using  a  front-type  matrix  inversion  that  assembles  a  portion  of  the 
matrix  and  solves  that  portion  before  assembling  the  next  portion  of  the  ma¬ 
trix.  The  front  solver's  efficiency  is  largely  independent  of  bandwidth  and 
thus  does  not  require  as  much  care  in  formation  of  the  computational  mesh  as 
do  traditional  solvers.  A  detailed  description  of  the  model  is  given  by 

Me Anally  et  al.  (1983). 

The  sediment  transp_ort_  mode_l_,  STUDH 

18.  Convection-diffusion  equat  ion.  The  sediment  transport  model,  STUDH 
solves  the  depth-integrated  convection-diffusion  equation  in  two  horizontal 
dimensions  fur  a  single  sediment  constituent.  The  basic  convection-diffusion 
equ.it  inn  is  presented  in  Ariathurai,  MacArthur,  and  Krone  (1977).  The  form  of 
the  solved  equation  is 


Where 

c  -  concentration,  kg/m* 

L  '  t  i me ,  sec 

u  flow  velocity  in  x-direetion,  m/sec 
^  primary  flow  direction,  m 


Figure  3.  Physical  model  limits 


proposed  30-ft  channels  leading  into  Baltimore  and  t  he*  existing  channels  else¬ 
where.  Channels  in  the  James  and  Elizabeth  Rivers  and  Thimble  Shoal  area  of 
the  Lower  Bay  were  molded  to  correspond  to  prototype  Information  collected  as 
late  as  19S1.  The  molded  area  of  the  model  extends  from  approximately  30  miles 
fft'shore  in  the  Atlantic  Ocean  to  the  heads  of  tide  for  all  tributaries  empty¬ 
ing,  into  the  Chesapeake.  The  entire  length  of  the  Chesapeake  and  Delaware 
:  <'&n)  Canal  and  a  portion  of  Delaware  Bay  are  also  molded.  Overhank  geometry 
is  reproduced  to  the  +20  ft  contour.  Model  limits  are  shown  in  Figure  3. 

12.  The  physical  model  was  designed  based  on  the  equality  of  Froude 
;  i he  r  s  ,  me, 1»»1  to  prototype,  reflecting  similitude  of  gravitational  effects. 
'Vm'  trio  scales  of  the  model  are  1:1,000  horizontally  and  1:100  vertically; 
t  he  Ji-t  rtion  ratio  is  10:1.  These  dimensions  and  Froudiar  model  laws  de- 
d  ’  following,  model -to-proto type  ratios: 


Character i st  ic 
Vertical  length 
Horizontal  length 
S  1  <>  pe 
2  i  m  e 

Ve  1  ni'  i  ty 
Vo! ume 
Di sc  ha rge 


Ra  t  i  o 
1:100 
1:1,000 
10:1 
1:100 
1  : 10 

1 : 100,000,000 
1 : 1,000,000 


ro-t  t  ype  ratio  for  salinity  is  1:1.  This  is  the  general  prac- 
r  !  ed -sc  1 1  e  mode  1  s  . 

model  was  designed  and  equipped  so  that  selected  prototype 
•  :  on>.  could  be  simulated  and  the  model  response  to  these  condi- 
.  \  detailed  discussion  ol  appurtenances  necessary  to  generate 

'  boundary  conditions  is  provided  in  Report  1  of  this  series 
"  14B3). 

Ilic  Mini ■  r  i c  a  1  Mod e  1  s 

i  I-  i-'!  a im  i  i  ..mdel  ,  RMA-2V ,  was  used  in  both  the 
i  :  ■  i  a  d  :  r  I  i  ■  t  b  t’iver  itudv.  The  noncohcs  i  ve  or  sand 


'  •  1 1 :  s por  t 


del,  STI'DH,  was  used  in  the  Thimble 


PART  !!: 


'Ml'  *'n|!i'!S 


Flio  ICvbrtd  vodelir\e  Techniijue 

■  1  - :  ;  I  1  •  -  !  inn  f  i  h  i  i  o  i  n  ■  i  t  s  i  A  in  t  ! .  i  s  s  l  " .  i  v  .  >: M '  !  r ^  <  1  .  ■■  :  .  r  \- 

'  I"  ’  tt*"  >-•«•  v:-iral  »»!*.•  1  with  a  ppr.  p  r  i  at  e  numerical 

Mi'  a:’:  i  -it  of  rat  ed  :i  I  1  •  d  t  r  pred  ir  t  i  channel  sed  i  :'i « ■  :i  t  at  ion  . 

:■  iMr.it  i  .mi  predictions  aro  cons;  iderod  to  he  more  r<*I  fable 
s  that  could  N-  made  using  either  physical  nr  numerical  models 

numerical  -ode  1  s  used  were  a  hvd  rod  yn  nr.  !c  mode  1  ,  called  RMA-?V  , 

*  r.t:v:p  .r t  model  ,  cal  1  ed  STl’DH.  The  Chesapeake  hay  physical 

the  nect  «sarv  water-surface  elevation  and  current  dat  a  to  drive 
vsir.i!  model  information  included  both  water-surf  ace  elevations 
all  t  '  "A-2V  boundaries  and  at  Locations  within  ea<-h  RVA-2V 

I  ■  1 1  and  Klicaheth  River).  These  locations  within  each  model 

d  in  i  i t  her  of  two  ways.  One  approach  is  to  use  these  loca- 

I I  hound -tries,  i.e.,  to  drive  each  RMA-2V  model  just  as  physi- 

■' in  drives  the  external  water  boundaries.  In  this  way,  the 
■nii-ic  model  is  real  1  V  used  to  interpolate  the  physical  model 

ludi  ill  A-2V  node  locations.  Another  approach  is  to  drive 

•  1  ■  ■  external  water  boundaries  and  use  the  internal  physical 

n  tor  comparison  with  RN  \-2V  results,  thus  nl  lowing  for  a  t  ra¬ 
nt  i  - '  i '  and  verification  of  t  ho  numerical  hydrodynamic  model, 
r*  1  ’  :  ve  sna  r  senes.s  of  internal  physical  model  data,  the  latter 

-  "  tor  this  studv  .  Ta'  resulting  currents  at  model  locations 
t  ’  used  hv  the  u  *- r  i  i  a  1  sediment  transport  model  ,  STl’Ml. 

lit  o’  t  h  ■  hvhrid  mode  1  i  ng  .approach  is  given  hv  Me  Anal  1  y  ct 

Cu.  j  Ipe.lke  1  *  v  I’hvsit  oil  'lode! 


i  . '  i  1  nl  flu-;  apeak-  i  •  1  oca!  f  d  on  Kent  Island  in 

c  ••  !•  I  i  •  o'  v  ;  ■  y,  d-hed  model  molded  in  i  nn- 

t  ’!  •-  r” r i ■  i  is  . u-iil  ‘.  I  ion.il  t'cean  Survey  (NOS)  charts.  At 

t  o'-.,  'll  "a  u .  r  ■  1  ■  1  p  i  1  si  ti  m  1  ••  1 1  ail  befit  mo  1  ded  with  the 


navigation  channel  using  existing  hybrid  modeling  techniques.  The  study  was 
designed  specifically  to  address  the  relative  changes  of  both  overall  shoaling 
and  distribution  of  shoaled  material  within  the  project  limits  caused  by 
deepening . 


8.  The  numerical  modeling  portion  of  the  hybrid  modeling  effort  con¬ 
sisted  of  two  separate  numerical  models,  which  were  referred  to  as  the  Thimble 
Shoal  model  and  the  Elizabeth  River  model.  Two  separate  models  were  needed  to 
encompass  the  study  area  rather  than  one  overall  model  because  of  the  varying 
nature  of  shoaled  material  along  the  project  navigation  channel.  The  sediment 
along  the  Thimble  Shoal  portion  of  the  navigation  channel  consists  predomi¬ 
nant!'''  of  noncohesive  material  while  the  sediment  along  the  Elizabeth  River 
portion  of  the  navigation  project  consists  primarily  of  clays  and  silts  (co¬ 
hesive  material).  Thus  the  noncohesive  version  of  the  numerical  sediment  model 
was  used  for  the  Thimble  Shoal  study  and  the  cohesive  version  of  the  numerical 
sediment  model  was  used  for  the  Elizabeth  River  study.  Sedimentation  in  a 
third  portion  of  the  overall  project,  referred  to  as  the  Atlantic  Ocean  Chan¬ 
nel,  is  not  addressed  in  the  main  body  of  this  report  but  is  presented  in 
Appendix  A. 


9 


Figure  9.  Elizabeth  River  water-surface  elevation  comparison  (Base  4) 


Figure  10.  Elizabeth  River  velocity  comparison  (Base  4) 


ihh!i'  for  each  t  imc— step  over  a  tidal  cu  If  wort*  then  stored  in  a  form  that. 

■  •on  Id  be  used  h  v  t  lie  sediment  mode)  . 

'id.  The  proposed  deepened  channel  was  simulated  by  deepening  the  rhannel 
irise  to  the  project  depth.  (Including,  an  additional  1  ft  for  maintenance  and 
dredging  tolerance)  without  changing  any  of  the  numerical  model  parameters. 

!  he  deepened  channel  model  was  then  run  using  the  four  sets  of  boundary  condi¬ 
tions.  Results  of  representative  station  velocity  comparisons  between  the 
ini  or  Leal  modal  and  the  physical  model  are  shown  in  Figure  11.  The  complete 
sot.  of  water-surface  elevation  and  velocitv  comparisons  for  al  1  test  condi¬ 
tions  are  shown  in  Plates  1 -?4.  The  hydrodynamic  model  plan  results  were 
store'*  i  :i  the  same  manner  as  the  base  results  for  use  in  the  sedimentation 


s.’d  i me r, tat  ion  simulation 

1  i .  in  order  to  properly  ad  just  the  Elizabeth  River  sediment  transport 
node! ,  representative  navigation  channel  shoaling  rates  had  to  be  established 
fro  iviilnhie  prototype  data.  To  that  end,  hydrographic  surveys  from  1 9f» 2  — 
were  acquired  from  the  I’S  Army  Engineer  District,  Norfolk.  The  method 
tsoi!  to  .•<>:- ptitc-  nnvftt.it  ion  channel  shoaling,  rates  was  to  compare  the  post- 
'  fudge  survey  with  the  following  prodr  edge  survey,  thus  establishing  the  amount 
>:  -inter  in  1  that  had  shoaled  during,  that  period.  For  the  existing,  4  t  -deep 
hansel  reaches  f Norfolk  "arbor  reach,  Cranev  Island  reach,  and  Lamberts  Fiend 
roach;  ,  live  such  o  spa r i sons  could  h,>  made ,  all  occurring  during  the  time 
•  :  i  •  id  *'r  »r  duly  l"7b  to  December  1981.  For  the  existing  40-1  t  -deep  channel 
reach  *  Dirt  No  r '■  > !  k  reach  ,  Town  Point  reach,  and  t  he  lower  and  Middle  readies 
'  the  Font,  "orn  dried),  five  comparisons  could  also  ho  made,  all  occurring 
vith;-'  t  h“  tii  e  period  >  r  >m  \ugust  i  97 1  to  February  1 980.  For  the  Vv-f  t-deep 
hanne  I  r>  -  h  )i  ( t  bo  Ppper  roach  of  t  he  Southern  branch),  six  comparisons  could 
'e  eh',  all  '  ic  •  ii  r  r  i  nr  within  the  period  from  Va  v  1  9b  1  to  Pebrunr  v  1  98 1  . 

i'.  flu  ti-i's  used  in  establishing  shoaling  rates  arc  as  follows: 


flh. a  nne  I 


■—ft  pro  iec  t 


Period 

df  Apr  i"Sl  to 
)  Dec  1  97"  t  e 


Channel 


Period 


40-ft  project  (1)  Nov  1975  to  Feb  1980 

(2)  Jul  1974  to  Sep  1975 

(3)  Feb  1974  to  Jul  1974 

(4)  Mar  1973  to  Sep  1973 

(5)  Aug  1971  to  Feb  1973 

35-ft  project  (1)  Mar  1979  to  Feb  1981 

(2)  Apr  1978  to  Mar  1979 

(3)  Feb  1974  to  Dec  1977 

(4)  Jun  1972  to  Sep  1973 

(5)  Oct  1969  to  Jun  1972 

(6)  May  1963  to  Oct  1969 

36.  To  allow  for  comparison  of  prototype  shoaling  with  the  sediment 
transport  model  shoaling,  the  existing  45-ft-deep  project  was  sectioned  as 
shown  in  Figure  12;  the  40-ft  project  was  sectioned  as  shown  in  Figure  13;  and 
the  35-ft  project  was  sectioned  as  shown  in  Figure  14. 

37.  The  prototype  shoaling  that  occurred  in  each  section  for  the  peri¬ 
ods  evaluated  and  the  averages  used  in  model  calibration  are  given  in  Table  1. 

38.  It  is  not  yet  economically  feasible  to  run  the  hydrodynamic  and 
sediment  transport  models  for  a  continuous  simulation  of  events  for  a  year  or 
even  several  months.  Therefore,  in  order  for  the  sediment  transport  model  to 
simulate  a  long  shoaling  period,  first  several  separate  events  were  simulated 
independently.  Then  each  event  was  extrapolated  over  the  appropriate  time  in¬ 
terval  and  all  were  combined  to  simulate  a  representative  annual  sedimentation 
cycle.  Four  separate  conditions  were  used  for  this  type  of  simulation  for  the 
Norfolk  Harbor  study  (identical  conditions  for  both  the  Elizabeth  Diver  model 
and  Thimble  Shoal  model) .  A  description  of  these  four  events  or  conditions  is 
given  in  paragraph  31.  The  impact  of  each  condition  on  the  long-term  shoaling 
in  the  model  was  determined  by  subdividing  the  year  of  events  into  the  por¬ 
tions  of  the  time  that  could  best  he  represented  by  a  particular  event. 

39.  The  combination  of  two  freshwater  inflows  and  two  tidal  ranges 
yields  four  events  or  conditions.  The  tidal  ranges  at  Old  Point  Comfort, 
which  corresponded  to  the  tidal  range  of  the  simulated  events,  were  2.40  and 
3.70  ft,  termed  the  mean  and  the  high  range  events.  The  proportion  of  the 
simulated  year  in  which  the  range  was  greater  than  3.05  ft  was  considered  the 
high  range  impact  period;  conversely,  the  time  period  proportion  less  than 
3.05  ft  was  the  mean. 

40.  These  duration  periods  were  ascertained  from  the  tidal  ranges 
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ype  shoaling  sections 


Figure  13.  Elizabeth  River  40-ft  project-prototype  shoaling  sections 


1 


Figure  1A.  Flizabeth 

River  35-ft  project-prototype  shoaling  sections 

generated  by  the  nine  major  astronomical  tidal  components  at  Old  Point  Comfort 

(Richards  and  Gulbrandsen  1982): 

Constituent 

Amplitude,  ft 

Period,  hr 

M2 

1.188 

12.A21 

N2 

0.265 

12.658 

S2 

0.230 

12.000 

K1 

0.170 

23.93A 

01 

0.1A6 

25.819 

V2 

0.051 

12.626 

PI 

0.0A8 

2A.066 

L2 

0.033 

12.192 

K2 

0.059 

11.967 

Since  actual  tides  were  not 

used,  meteorological 

effects  were  neglected; 

therefore  the  total  energy 

may  be  slightly  lower 

than  actual  values.  At  Old 

Point  Comfort,  the  average  astronomical  tidal  range  was  2.AA  ft;  the  mean 
simulation  event  had  a  range  of  2. AO  ft.  The  high  range  events  would  then  be 


simulated  16  percent  of  the  model  year  and  the  mean  range  events,  8A  percent. 
41.  The  Inflow  conditions  in  the  physical  model  were: 


Mean  Flow,  cfs 


Nansemond  River 

372 

1,063 

Chickahominy  River 

289 

826 

Appomattox  River 

967 

2,763 

James  River 

7,249 

20,711 

Total  inflow  to 

70,000 

200,000 

Chesapeake  Bay 


42.  The  duration  of  flows  in  the  James  River  near  Richmond  are  from  the 
Virginia  Department  of  Water  Resources  (1974).  From  this  reference,  it  was 
calculated  that  the  high  inflow  event  (greater  than  13,980  cfs  in  the  James 
River  near  Richmond)  was  exceeded  only  10  percent  of  the  time.  Accordingly, 
the  high  flow  events  would  be  modeled  for  10  percent  of  the  time  and  the  mean 
flow  events  the  remaining  90  percent. 

43.  The  combination  of  these  independent  inflow  and  range  events  would 

be  : 


Duration 


Event 

Range 

Inflow 

Percent  of  Total 

1 

High 

High 

1.6 

2 

Mean 

High 

8.4 

3 

High 

Mean 

14.4 

4 

Mean 

Mean 

75.6 

44.  The  various  model  parameters  established  during  calibration  were: 


Time-step  length 

Crank  Nicholsen  implicitness 
factor 

Manning’s  n  roughness 

Critical  shear  stress  for 
deposition 

Critical  shear  stress  for 
erosion 
Top  layer 
Second  layer 

Particle  erosion  constant 

Effective  settling  velocity 
James  River 

Within  F.lizabeth  River 
Dry  weight  density  of  deposit 

( Continued) 


1,800  sec  (0.5  hr) 

0.70 

0.017 

0.02  N/m2 

0.02  N/m2 
0.06  N/m2 

0.0012  kg/m2/sec 

0.0008  m/sec 
0.0001  m/sec 

500  kg/m^ 
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Sediment  Concentration  Boundary  Conditions* 


High  Freshwater 

Mean  Freshwater 

Inflow 

Inflow 

James  River  (upstream) 

0.050  kg/m3 

0.035  kg/m3 

Chesapeake  Bay 

'i 

(downstream) 

0.010  kg/nr 

0.010  kg/nK 

Upper  Elizabeth  River 

0.001  kg/m3 

0.001  kg/ra 

*  Values  based  on  Onishi  and  Wise  (1978). 


45.  The  length  of  the  time-step  was  chosen  to  not  allow  sediment  to 
pass  completely  through  an  element  during  one  time-step.  Many  of  the  other 
parameters  are  present  to  address  the  cohesive  sediment  transport  processes  of 
this  region.  Sediment  within  the  Elizabeth  River  was  largely  silts  and  clays 
and  so  was  cohesive  in  character.  The  critical  shear  stresses  for  deposition 
and  erosion  were  chosen  to  allow  deposition  and  resuspension  (erosion)  in  the 
Norfolk  Harbor  reach  and  primarily  deposition  alone  throughout  most  of  the 
Elizabeth  River. 

46.  A  model-to-prototype  comparison,  model  calibration,  is  shown  in 
Figure  15.  These  model  shoaling  rates  were  the  superposition  of  the  influence 
of  each  of  the  four  flow  events. 


Figure  15.  Elizabeth  River  model  sedimentation  verification 
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47.  The  model's  maximum  shoaling,  about  270  mm/ yr  (prototype  maximum  of 
480  mm./vr)  occurred  in  the  sane  area  of  the  Norfolk  Harbor  reach  as  in  the 
prototype.  The  model  and  prototype  shoaling  rates  dropped  within  the  Elizabeth 
River  channels.  The  Craney  Island  disposal  area  effectively  blocked  material 
from  entering  the  Elizabeth  River  mouth  from  the  west;  also  the  currents  are 
somewhat  stronger  due  to  the  smaller  cross  section.  Both  prototype  and  model 
results  indicated  a  minimum  deposition  around  Lamberts  Bend. 

4k.  In  tin  Port  Norfolk  reach,  the  model  showed  a  peak  shoaling  rate  of 
about  43  m:n.'vr ;  this  shoaling  rate  showed  a  decline  from  there  to  the  upper 
to  .!•.'(  of  the  project.  The  prototype  data  showed  shoaling  rates  in  Elizabeth 
m.o  r  wi greater  variation  than  those  in  the  model,  and  with  a  peak  in  the 
;  min  basin  at  the  head  of  the  upper  reach.  The  model  sediment  source  was 
'1:1.0.  :  1  /  fro,:,  t  h<  James  River,  while  the  prototype  might  have  had  a  local 

in  t.:!  it  ion.  However,  the  overall  model  shoaling  pattern  was  quite 
’  •  i-  it  of  prototype  and  the  magnitudes  were  quite  reasonable. 

Thimble  Shoal 


o  !  r  1  i  •  si  mu  1  a  t  i  on 

■>.  .he  sane  procedures  used  in  the  Elizabeth  River  simulations  were 
.0  :  i  ?  Thimble  Shoal  simulations.  The  boundary  file  source  was  still  the 

o'  i  ,■  i  1  "ode  1  data  and  .all  station  comparisons  were  made  to  the  physical  model 
i  •  a  .is  in  t.  tie  Elizabeth  River  simulations. 

'••').  \'os  chart  numbers  12248  and  1  2222  were  used  as  the  basis  for  the 
••'■■sis  overlay  representation  of  the  lower  portion  of  Chesapeake  Bay  and  the 
on  1  r.an-e  to  the  .Tames  River  estuary.  The  Thimble  Shoal  mesh  constructed  con¬ 
sisted  of  1,393  nodes  and  428  elements.  Figure  lb  shows  the  resulting  mesh. 

The  ups  t.  ream  and  downstream  water  boundary  conditions  of  the  Thimble  Shoal 
mesh  were  velocity  or  discharge  at  Hampton  Roads,  the  entrance  to  the  Elizabeth 
Ri ver ,  and  the  northern  cutoff  of  the  Chesapeake  Bay  and  tidal  ole vat  ions  at 
*  he  ocean. 

hi.  Iho  sane  type  of  physical  model  data  furnished  for  the  Elizabeth 
River  simulations  was  also  used  for  the  Thimble  Shoal  simulations.  All  the 
'■■st  conditions  remained  the  same,  on  1 v  the  computational  mesh  and  boundary 
locations  were  different.  Figures  17  and  18  show  the  locations  of  the  c ompa ra¬ 
ti  y.  physical  model  stations  with  numerical  model  nodes. 
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Figure  31.  Thimble  Shoal,  sedimentation  results, 
annual  rate  comparison 

ion  between  ranges  -3b, 000  to  -40,000  ft  was  reduced  from  about  90  to 
m/yr.  Tt  is  possible  that  the  lack  of  this  material  source  is  responsibl 
tin.  apparent  shoaling  reduction  in  this  region. 

7P.  The  amount  of  shoaling  represented  an  annual  volume  increase  of 
t  ’!1  percent  from  base  to  deepened  plan  conditions.  If  one  considered  tin 
sedimentation  rates  (scouring  and  shoaling),  the  plan  demonstrated  a 
•  rcent  i no  reuse  over  base,  the  first  determination  (20  percent)  would  be 
appropriate  manner  to  calculate  the  potential  increase  maintenance  dredg- 
volume.  as  areas  that  showed  net  scour  would  not  need  to  be  dredged. 


Figure  10.  Elizabeth  River,  sedimentation  results,  annual  rate  comparison 

deposition,  each  of  the  four  conditions  for  Thimble  Shoal  consisted  of  both 
scour  and  deposition.  The  conditions  are  related  in  that  deposition  in  one 
condition  is  erosion  in  another  condition,  making  results  of  a  single  condi¬ 
tion  less  meaningful  than  in  the  Elizabeth  River  simulation.  For  this  reason 
only,  the  yearly  base  versus  plan  comparison  Is  shown  (Figure  30)  for  Thimble 
Shoal  .  The  channel  deepening  had  its  greatest  impact  on  the  lower  four  reaches 
(Norfolk  Harbor,  Craney  Island,  Lamberts  Bend,  and  Port  Norfolk  reaches). 
Velocities  in  this  region  were  large  enough  to  resuspend  some  of  the  freshly 
deposited  material;  therefore  the  deepenings  reduction  of  velocity  here  had  a 
greater  effect.  The  annual  volume  of  channel  deposition  was  23  percent  higher 
over  the  length  of  the  channel  for  plan  relative  to  base. 

Thimble  Shoal  Channel 


68.  Comparison  of  model  results  between  base  and  after  deepening  (plan) 
conditions  is  shown  in  Figure  31.  The  shoaling  peaks,  both  west  and  east  of 

f  ho  Bridge,  moved  outward  toward  the  ends  of  the  channel.  This  would  seem  to 
indicate  that  the  material  depositing  entered  from  each  end,  and  with  the 
deepened  conditions  the  current  velocities  were  not  capable  of  carrying  the 
material  as  far  into  the  channel. 

69.  The  seaward  peak  was  increased  from  about  30  to  70  mm/yr.  The  peak 
landward  shoaling  rate  was  reduced  from  about  63  to  43  mm/yr.  Scour  in  the 
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Elizabeth  River,  sedimentation  results,  base  to  plan  (Flow  3) 
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rigure  26.  Elizabeth  River,  sedimentation  results,  base  to  plan  (Flow  1) 
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Figure  27.  Elizabeth  River,  sedimentation  results,  base  to  plan  (Flow  2) 


PART  TV: 


S  F  D  T M  E  NT  A  T  T  0  N  RESENTS 


Elizabeth  River 

64.  Comparisons  of  existing  (hase)  and  deepened  (plan)  channel  sedi¬ 
mentation  results  for  each  flow  event  are  shown  in  Figures  26-79.  These  fig¬ 
ures  were  derived  on  an  annual  basis  for  each  event,  so  the  magnitudes  have 
only  relative  significance.  Fach  figure  also  includes  the  average  suspended 
sediment  concentration  profile  along  the  channel.  The  suspended  sediment  con¬ 
centration  profile  indicated  a  slight  reduction  in  the  channel  for  the  plan 
for  all  four  events.  There  is  an  increase  in  shoaling  rate  in  most  reaches  of 
the  Flizabeth  River  under  all  four  conditions  of  the  new  deepened  channel. 

65.  The  rate  of  sedimentation  was  highest  in  the  Norfolk  Harbor  reach 
for  all  four  events  for  both  plan  and  base  conditions.  However,  the  high 
range  events  (Events  1  and  1)  shoaled  substantially  less  in  this  reach  than 
did  the  lower  range  events  (Events  2  and  4).  The  greater  tidal  prism  of  the 
Flizabeth  River  produced  higher  flow  magnitudes  and  thus  was  able  to  reduce 
deposition.  The  difference  between  the  plan  and  base  shoaling  rates  was  great¬ 
est  for  the  high  range  events.  The  increased  channel  depth  reduced  current 
velocity  magnitudes  and  subsequently  the  bed  shear  stress  and  therefore  did 

not  allow  as  much  resuspension  of  freshly  deposited  material.  For  the  lower 
tidal  range  events,  the  bed  shear  stress  did  not  generate  a  great  deal  of  re¬ 
suspension  in  base  and  so  the  increased  channel  depth  did  not  increase  net 
deposition  very  much. 

66.  Another  interesting  observation  is  that  the  high  flow  events 
(Events  1  and  2)  shoaled  somewhat  less  than  the  lower  flow  event.  Tt  should 
be  noted  that  the  events  termed  "low"  for  tidal  range  and  freshwater  inflow 
art*  in  fact  "mean"  events.  Tt  would  be  a  mistake  to  place  a  great  deal  of  em¬ 
phasis  on  such  a  detailed  analysis  of  model  results  without  substantiation  from 
the  prototype.  Since  the  model  was  calibrated  to  a  large  time  period,  com¬ 
bining  a  mulfitudr  of  events,  its  validity  would  be  limited  to  that  degree  of 

i'  nm pa  r  i  son  . 

F7.  Figure  IF  shows  the  comparison  utilizing  a  combination  of  the  four 
events  for  an  annual  estimate  of  shoaling  for  base  and  plan  conditions  by 
,  Funnel  reaches.  The  base  versus  plan  comparisons  for  each  of  the  four  condi- 
t  ions  are  not  displayed.  Fn 1 i ke  the  Flizabeth  River  results  which  were  mainly 


this  particular  region.  The  purpose  of  the  study,  to  determine  the  maintenance 
dredging  change  due  to  channel  deepening,  should  still  be  addressable  since 
deepening  of  the  channel  should  not  result  in  the  channel  trapping  any  more 
material  from  the  migratory  shoal. 

62.  The  model  results  showed  the  same  shoaling  pattern  as  that  of 
prototype  but  were  generally  lower  in  magnitude. 

63.  Shoaling  outside  the  channel  domain  was  not  noticeable  in  the 
sediment  transport  model;  this  seems  to  be  in  line  with  a  previous 
estimate  of  3.7  mm/yr  in  this  region  (Ludwick  1981). 
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took  place  upstream  of  the  dredged  navigation  channel  around  the  narrows  at 
the  James  River  mouth. 

60.  The  model  calibration  results  from  the  combination  of  the  three 
simulated  events  are  shown  in  Figure  25.  The  largest  prototype  shoaling  rate 
occurred  between  ranges  -20,000  and  -24,000  ft  and  is  about  100  mm/yr.  The 
sedimentation  model  showed  a  peak  between  ranges  -24,000  and  -28,000  ft  and  was 
about  65  mm/yr.  The  model  shoaling  peak  here  is  due  to  the  ebb  plume  from  the 
James  River.  Another  interesting  feature  is  the  net  scour  indicated  by  the 
surveys  in  the  segment  from  range  -36,000  to  -40,000  ft  of  about  60  mm/yr. 

This  may  be  an  indication  that  the  system  is  still  adjusting  for  the  channel 
presence  or  perhaps  the  location  of  the  natural  channel  is  shifting.  The 
model  indicated  a  scour  of  about  90  mm/yr  at  this  segment. 


UPSTREAM 

DOWNSTREAM 

Figure  25.  Thimble  Shoal  sedimentation  verification 

61.  The  second  shoaling  peak  in  the  prototype  (between  ranges  +12,000 
and  +16,000  ft)  was  about  60  mm/yr.  This  shoal  was  in  the  near  vicinity  of 
the  "Tail  of  the  Horseshoe"  shoal,  a  mobile  shoal  just  north  of  the  channel 
which  according  to  the  findings  of  Ludwick  (1981)  has  "experienced  net  accre¬ 
tion  on  its  southern  flank  during  the  past  century."  This  could  be  an  ex¬ 
planation  for  the  second  peak.  The  model  showed  a  less  distinct  peak  than  the 
prototype;  this  occurred  between  ranges  +16,000  and  +20,000  ft  and  was  about 
30  mm/yr.  The  sedimentation  in  the  model  is  low  here  in  spite  of  using  an  ef¬ 
fective  sediment  diameter  for  transport  smaller  than  was  generally  found  in 
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Two  other  periods  of  shorter  duration  (less  than  1  year)  were  available  for 
analysis;  but  due  to  the  very  low  shoaling  rates  in  the  Thimble  Shoal  Channel 
and  the  limited  accuracy  of  hydrographic  surveys,  the  decision  was  made  to 
omit  these  short  duration  comparisons  and  use  only  the  1975-1981  period  for 
calibration  of  the  model. 

55.  To  allow  for  comparison  of  prototype  shoaling  with  the  sediment 
transport  model  shoaling,  the  Thimble  Shoal  Channel  was  sectioned  as  shown  in 
Figure  24. 


THIMBLE  SHOAL  CHANNEL 

Figure  24.  Thimble  Shoal  Channel-prototype  shoaling  sections 

56.  The  prototype  shoaling  that  occurred  in  each  section  for  the  period 
evaluated  is  given  in  Table  2. 

57.  The  simulation  of  an  annual  event  was  conducted  in  the  same  manner 
as  described  in  paragraph  38.  However,  difficulties  arose  in  the  hydrodynamic 
simulation  of  F.vent  3.  Therefore  Fvent  1  was  used  for  both  high  range  events. 

58.  The  length  of  the  sedimentation  model  time-step  selected  was 
1,800  sec  (0.5  hr).  This  was  small  enough  not  to  allow  advection  of  sediment 
completely  through  an  element  during  one  time-step.  The  computational  grid  was 
identical  with  that  of  the  hydrodynamic  model.  Other  model  parameters  were 
generally  as  follows: 


Crank  Nicholsen  implicitness  factor 
Manning's  n  roughness 

effective  particle  diameter  for  transport 
effective  settling  velocity 
Boundary  concentration 


0.65 
0.017 
0.05  mm 
0.0008  m/sec 
0.010  kg/m 


59.  Tii  i  s  modeling  effort  was  conducted  assuming  noncohesive  sediment 
transport  with  the  given  effective  particle  diameter  for  transport  as  being 
-nurse  silt  to  very  fine  sand.  This  transport  size  appears  most  applicable  in 
the  western  portion  of  the  navigation  channel  (refer  to  Figure  4  from  Ludwick 
M  9H 1  ) ) .  Fven  using  this  fairly  small  diameter,  most  of  the  sediment  movement 
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Figure  23.  Thimble  Shoal  numerical  model-physical  model 
velocity  comparison  (Plan  1) 

model  and  the  physical  model  are  shown  in  Figure  23.  Water-surface  elevation 
and  velocity  comparisons  for  all  test  conditions  can  be  found  in  Plates  25-48. 
After  verification  was  achieved,  the  hydrodynamic  model  results  were  stored 
for  use  in  the  sedimentation  simulations. 

Sed iment  simulation 

54,  In  order  to  properly  adjust  the  Thimble  Shoal  sediment  transport 
model,  representative  navigation  channel  shoaling  rates  had  to  be  established 
from  available  prototype  data.  Hydrographic  surveys  from  1960-1981  were 
acquired  from  the  Norfolk  District.  The  method  used  to  compute  navigation 
channel  shoaling  was  the  same  as  that  described  in  paragraph  34  for  the 
Fllzaheth  River  model.  The  single  period  from  10  February  1975  to  14  April 
1 9H 1  served  as  the  calibration  standard  for  the  sediment  transport  model. 


Figure  21.  Thimble  Shoal  water-surface  elevation  comparison  (Base  4) 
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V/.  (’a  ]  i  hr  at  i  on  of  the  numerical  hydrodynamic  simulation  proceeded  in 

t  he  same  fashion  as  that  of  the  Kl izabeth  River  simulation.  The  friction  or 
odd v  ili ftusivity  parameters  were  ad  justed  until  satisfactory  agreement  was 
achieved  between  the  numerical  and  physical  models.  Then  the  numerical  model 
was  verified  against  Base  1  conditions.  Finally,  Base  2  and  3  conditions  were 
run  to  complete  the  series.  Figures  19-22  show  tidal  and  velocity  comparisons 
for  critical  representative  locations.  The  hydrodynamic  results  were  then 
stored  for  use  in  the  sedimentation  simulations. 

93.  The  nodes  representing  the  Thimble  Shoal  Channel  were  deepened  to  a 
uniform  depth  of  S8  ft.  This  plan  depth  included  the  project  depth  and  3  ft 
for  dredging  tolerance  and  advance  maintenance.  The  same  boundary  conditions 
used  for  t  lie  base  simulations  were  then  used  for  the  deepened  channel  config¬ 
uration.  Results  of  a  representative  station  comparison  between  the  numerical 
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in  i  Mill  i  '  mOmi, 


figure  19.  Thimble  Shoal  water-surface  elevation  comparison  (Base  1) 


PART  V:  CONCLUSIONS 


71.  Based  on  sedimentation  results  from  the  Elizabeth  River  numerical 
model,  the  increase  in  annual  shoaling  caused  by  channel  deepening  as  proposed 
will  be  21  percent.  The  distribution  of  shoaled  material  will  not  be  signifi¬ 
cantly  altered,  cither  than  a  slight  increase  in  skewness  toward  the  downstream 

end  . 

72.  Based  on  sedimentation  results  from  the  Thimble  Shoal  numerical 
model,  the  increase  in  shoaling  caused  by  channel  deepening  as  proposed  will 
be  about  20  percent.  The  distribution  of  shoaled  material  will  be  slightly 
altered  in  that  both  the  upper  and  lower  channel  shoaling  peaks  which  presently 
exist  will  tend  to  migrate  even  more  toward  the  ends  of  the  dredged  channel. 

71.  The  estimate  of  shoaling  for  the  new  Atlantic  Ocean  Channel,  de¬ 
scribed  in  the  appendix  to  this  report,  is  about  200,000  cu  yd  annually. 
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APPENDIX  A:  ATLANTIC  OCEAN  CHANNEL  SEDIMENTATION  PREDICTIONS 


Introduction 


ground 


1,  The  development  of  the  Atlantic  Ocean  Channel  is  part  of  the 
iosed  improvements  to  channels  approaching  Norfolk  Harbor  and  anchorages, 
i  new  channel,  located  off  Virginia  Beach,  is  to  be  57  ft  below  mean 
water  and  have  a  width  of  1,000  ft  over  a  length  of  10.6  statute  miles, 
s  stated  depth  of  the  channel  is  the  project  depth.  The  actual  depth 
the  proposed  new  channel  with  the  dredging  tolerance  and  advance  main- 
ince  will  be  60  ft  below  mean  low  water.  Figure  A1  shows  the  project 
nils . 

- -  37° 

CHESAPEAKE  BAY  ENTRANCE 


-CAPE 

HENRY 
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Figure  Al.  Project  location 


2.  The  purpose  of  this  portion  of  the  broader  study  was  to  analyti- 
ly  or  empirically  predict  expected  shoaling  rates  in  this  channel  once  it 
developed.  This  was  to  be  accomplished  using  data  that  were  available. 


Al 


e  greater  than  30  cm/sec  9.30  percent  of  the  time.  For  the  60-ft  Atlantic 
ean  Channel  the  wave  orbital  velocities  are  greater  than  30  cm/sec  8.39  per- 
nt  of  the  time.  Therefore  the  ability  of  the  waves  to  at  least  place  fine 
nd  into  suspension  would  only  be  changed  by  approximately  1  percent  with 
e  deepened  channel. 
drographic  survey  comparisons 

32.  At  the  site  of  the  proposed  Atlantic  Ocean  Channel  three  hydro- 
aphic  surveys  can  be  compared  to  investigate  gross  changes  in  bathymetry. 
ie  1947  survey  is  a  patchwork  with  many  areas  of  no  data.  The  1980  survey 
tcluded  the  entire  study  area.  The  1969  survey  did  not  extend  past  75°49' 
.ongitude)  and  36°31 '  (latitude);  therefore  only  the  first  7.5  miles  of  the 
lannel  can  be  compared. 

33.  Both  the  1969  and  1980  survey  data  were  collected  by  an  on-line 
■quisition  system  (Gardner  1982*)  and  referenced  to  mlw. 

34.  To  compare  surveys,  comparison  points  were  located  on  the  hydro- 
raphic  surveys  and  the  10  depths  surrounding  this  point  were  recorded  and 
^eraged.  From  the  23  comparisons  made  between  the  1969  to  1980  data  there 
as  a  natural  deepening  or  erosion  that  averaged  about  20  mm/ yr .  From  the 
947  to  1969  surveys  the  computed  erosion  averaged  about  120  mm/yr. 

35.  The  1969  and  1980 
urveys  were  compared  in  more 
etail  using  a  computer-based 
ystem  developed  at  the  US 
rmy  Engineer  Waterways  Ex- 
eriment  Station  (LaGarde  and 
eltzel  1980).  Figure  A9 
hows  the  coverage  of  the 
urveys  compared.  The  1969 
urvey  sheets  were  reduced  to 
he  scale  of  the  1980  sheets, 
oth  surveys  were  contoured  by 
and  and  then  digitized  and 

Figure  A9 .  Hydrographic  survey  coverage 


_  • 

Personal  communication,  LT  J.  C.  Gardner,  Jr.  (1983),  N0AA,  Norfolk,  Va . 
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Figure  A6 .  Chesapeake  Bay  model  data  (Test  2) 
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NORFOLK  HARBOR  AND  CHANNELS 
CHESAPEAKE  BAY  MODEL  DATA 
AC0002  (TEST  1) 


Condition  2:  The  total  bay  discharge  remained  at  200,000  cfs;  how¬ 
ever,  the  repetitive  cosine  tide  range  was  changed  to  +1.50  ft, 
which  approximates  a  neap  tide  condition. 

Condition  3:  The  total  bay  discharge  was  70,000  cfs,  which  simu¬ 
lates  the  long-term  average  flow  into  the  Chesapeake  Bay  from 
all  its  tributaries.  The  approximate  spring  tide  range  of 
+2.4  ft  was  generated  at  the  ocean. 

Condition  4:  The  total  bay  discharge  remained  at  70,000  cfs  and 

the  approximate  neap  tide  range  of  +1.5  ft  was  generated  at  the 
ocean . 

28.  Velocity  plots  (Figures  A5-A8)  of  the  data  indicate  flood  pre¬ 
dominance  at  all  depths  except  the  surface  and  this  pattern  existed  during  all 
the  test  conditions  (flood  values  are  positive).  An  analysis  of  direction 
data  indicated  that  during  a  majority  of  the  four  test  conditions  the  flow  was 
along  the  axis  of  the  channel.  The  spring  tidal  range  conditions  produced  a 
few  velocities  greater  than  1  fps;  however,  during  the  neap  tidal  range,  all 
velocities  were  less  than  1  fps.  This  velocity  would  be  the  minimum  required 
to  move  sediment  found  in  the  study  area  from  the  bed  into  suspension.  Results 
from  the  test  with  the  Atlantic  Ocean  Channel  installed  indicated  the  same 
general  trends  but  the  velocities  were  reduced.  Therefore  tidal  currents  alone 
would  not  be  sufficient  to  cause  any  significant  change  in  the  bed. 

Wave  climatology  study  for  Chesapeake 
Bay  entrance  and  Virginia  Beach _ 

29.  Wave  climatology  statistics  of  the  Chesapeake  Bay  entrance  and 
Virginia  Beach  areas  for  the  20-year  period  from  1956-1975  are  given  in 
Tables  A1  and  A2,  respectively  (Jensen  1982*).  These  wave  statistics  (sig¬ 
nificant  waves  in  height  and  period)  include  20-year  summaries  (percent  oc¬ 
currence)  for  the  entire  period  for  all  directions. 

30.  From  the  20-year  summary  of  wave  statistics  for  Virginia  Beach,  it 
is  observed  that  from  58,440  significant  waves,  81  percent  had  a  height  of 
1.49  m  or  less  and  71.15  percent  of  the  periods  were  7.9  sec  or  less.  The 
largest  significant  wave  was  5.35  m  and  the  average  was  0.55  m. 

31.  Using  the  method  outlined  by  Beauchamp  (1974)  to  determine  the 
maximum  orbital  velocity,  it  was  determined  that  at  the  present  average  depth 
of  the  location  of  the  proposed  channel  (56.96  ft)  wave  orbital  velocities 


*  Personal  communication,  R.  E.  Jensen,  1982,  Wave  Information  data  base, 
provided  by  Coastal  Engineering  Research  Center,  USAEWES. 


dredged  channel,  located  at  Port  Said  Harbor,  Egypt.  In  the  method  outlined, 
an  average  steady  current  characterized  by  a  velocity  and  a  wave  height  and 
period  is  required.  From  this  information  a  maximum  oscillatory  current  at 
the  bed  due  to  wave  action  is  calculated.  Values  are  obtained  for  the  two 
components  of  sediment  transport  (bed  and  suspended  load),  based  on  the 
maximum  oscillatory  current  and  selected  descriptive  sediment  and  water 
parameters,  in  addition  to  other  laboratory  derived  variables.  In  estimating 
the  deposition  rate,  it  is  assumed  that  part  of  the  suspended  load  and  all 
of  the  bed  load  are  deposited  within  the  channel.  The  change  in  deposition  is 
based  on  the  changes  in  the  load-carrying  capacity.  Thus  the  deposition  is 
the  sum  of  the  bed  load  and  difference  between  the  suspended  load  approaching 
and  leaving  the  channel. 

26.  Ludwick  (1981)  proposed  a  short  method  that  produces  reasonable 
results.  He  simplifies  the  equation  for  deposition  of  fine-grained  sediment 
in  the  presence  of  a  moving  tidal  current  by  assuming  that  for  a  specified 
short  period  of  time,  the  shear  stress  at  the  bed  is  substantially  less  than 
the  critical  shear  stress  for  deposition.  Thus  the  mass  of  sediment  deposited 
per  unit  time  is  equal  to  the  average  concentration  times  an  average  settling 
velocity.  A  thickness  rate  of  deposition  is  determined  from  the  mass  of 
sediment  deposited  by  estimating  an  in-place  bulk  density  of  the  sediment.  The 
deposition  period  referred  to  is  taken  as  a  1-hr  period  at  slack  water  before 
ebb  and  at  slack  water  before  flood.  This  assumption  results  in  deposition 
only  occurring  2.0  hr  out  of  every  12.42  hr. 

Results 


Chesapeake  Bay  physical 
model  data  analysis 

27.  Velocity  measurements  were  made  in  the  Chesapeake  Bay  physical 
model  at  a  station  east  of  Rudee  Inlet  and  approximately  in  the  center  of  the 

proposed  Atlantic  Ocean  Channel.  This  station  is  designated  as  AC0002.  This 

was  done  to  examine  the  impact  of  channel  deepening  on  current  velocities. 

The  following  were  the  four  steady-state  boundary  conditions: 

Condition  1:  The  model  was  operated  under  a  repetitive  cosine  tide 
having  a  range  of  +2.40  ft  at  the  Atlantic  Ocean  control  sta¬ 
tion.  This  approximates  a  spring  tide  condition  in  Chesapeake 

Bay.  The  total  bay  freshwater  river  inflow  was  a  constant 
200,000  cfs,  which  represents  a  relatively  high  flooding 
condition. 


Available  analytical 
and  empirical  methods 

21.  Many  techniques  have  been  used  in  the  past  to  predict  the  effects 
of  channel  enlargement  on  channel  shoaling.  These  range  from  rule-of-thumb 
predictions  to  elaborate  physical  model  investigations.  During  the  past  few 
years,  there  has  been  progress  in  the  use  of  numerical  models  or  combined 
physical  and  numerical  models  (hybrid  studies)  for  making  these  predictions. 
Empirical  techniques  and  analytic  solutions  have  been  used  extensively  and 
provide  valuable  easily  obtained  solutions;  they  are,  in  most  cases,  scien¬ 
tifically  based. 

* 

22.  Trawle  has  done  a  thorough  review  of  some  of  the  available  analytic 
and  empirical  methods  for  estimating  infill  rates  in  offshore  channels.  These 
techniques  can  be  used  when  consideration  is  being  given  to  a  channel  enlarge¬ 
ment  or  to  the  use  of  advance  maintenance. 

* 

23.  Trawle  describes  four  analytically  based  and  three  empirically 
based  methods.  He  identifies  the  four  analytic  methods  as  the  Moriches  Inlet 
method,  the  Lean  method,  the  Oregon  Inlet  method,  and  the  Lamble  method.  The 
three  empirical  methods  discussed  are  identified  as  the  Gole  method,  the 
Simplified  Shoaling  Rate  method,  and  the  Volume  of  Cut  method.  These  repre¬ 
sent  some  of  the  better  documented  techniques  currently  being  used. 

24.  Boicourt  (1981)  uses  what  may  be  termed  a  single-grain  character¬ 
ization.  This  technique  is  implemented  by  collecting  near-bottom  velocity 
data  and  determining  a  mean  grain  size  for  the  sediment  in  the  area.  The  cur¬ 
rent  meter  data  3  m  or  so  above  the  bottom  are  assumed  to  be  appropriate 

for  a  conservative  estimate  of  the  threshold  velocity  for  initiation  of  grain 
movement.  A  threshold  velocity  for  initiation  of  grain  movement  for  that 
particular  sediment  size  is  calculated.  Then  he  determines  from  his  velocity 
record  how  much  of  the  time  this  value  is  exceeded.  By  comparing  these  data 
he  makes  a  qualitative  assessment  of  when  sediment  transport  events  are  likely 
to  occur.  This  comparison  could  be  done  seasonally. 

25.  Kadib  (1976)  describes  yet  another  method  to  predict  sedimentation 
at  offshore  dredged  channels.  His  approach  is  based  on  theoretical  studies 
and  experiences  gained  with  maintenance  dredging  at  the  Suez  Canal  offshore 


*  M.  J.  Trawle.  Engineer  Technical  Letter  1110-2-293,  Engineering  and  De¬ 
sign,  "Shoaling  Predictions  in  Offshore  Navigation  Channels,  Analytical  and 
Empirical  Methods.” 


produce  a  net  sediment  transport.  Actual  near-bottom  velocities  have  a 
larger  forward  velocity  of  shorter  duration  under  the  wave  crests  and  a 
smaller  backward  velocity  of  longer  duration  under  the  troughs.  A  net  sedi¬ 
ment  transport  can  be  caused  by  this  small  difference  between  large  quantities 
of  sediment  in  motion.  Laboratory  observations  and  field  work  have  confirmed 
that  sediment  transport  outside  the  breaker  zone  is  composed  of  both  bed 
and  suspended  load.  The  near-bottom  shear  velocity  is  the  same  order  of 
magnitude  as  the  fall  velocity  of  the  sediment  particles. 

19.  Hales  (1980)  discusses  in  detail  the  combined  effects  of  current 
and  wave  motion.  He  notes  that  waves  change  characteristics  as  they  progress 
from  relatively  calm  water  into  regions  of  streaming  water.  Waves  traveling 
with  the  current  experience  an  increase  in  length  and  celerity  and  a  decrease 
in  height.  Waves  traveling  against  the  current  increase  in  height  and  de¬ 
crease  in  length  and  celerity  until  a  limiting  steepness  occurs  that  depends 
both  on  the  initial  wave  characteristics  and  the  strength  of  the  opposing 
current  (Hales  1980). 


Review  of  Potential  Solution  Techniques 

Introduction 

20.  In  the  WESHE  proposal  dated  20  August  1981,  Subject:  Time  and  Cost 
Estimates  for  Norfolk  Harbor  55-Foot  Channel  Study,  a  plan  of  study  for  the 
Atlantic  Ocean  reach  was  proposed:  "The  Atlantic  Ocean  area  of  the  Chesapeake 
Bay  model  has  not  been  verified  to  properly  reproduce  prototype  currents; 
therefore  the  Atlantic  Ocean  Channel  is  not  amenable  to  a  hybrid  modeling 
technique  employing  the  bay  model.  It  is  possible  to  perform  a  sedimenta¬ 
tion  study  of  that  area  by  means  of  a  completely  numerical  modeling  effort, 
including  modeling  of  wind  waves  and  longshore  currents,  but  is  considered 
inadvisable  because  the  cost  (including  extensive  field  data  collection)  does 
not  appear  to  be  justified  by  the  magnitude  of  the  shoaling  problem  in  that 
channel.  Instead  we  recommend  an  analytic  study  using  available  data  to 
predict  an  expected  channel  shoaling  rate.  Presently,  we  have  an  active  re¬ 
search  study  under  way  to  develop  the  most  effective  procedures  to  address 
sediment  predictions  in  entrance  channels.  The  most  effective  technique 
resulting  from  this  research  effort  will  be  used  to  investigate  this  portion 
of  the  study." 


Theoretical  Background  of  the  Problem 


14.  A  thorough  treatment  of  the  problems  of  unsteady  oscillatory  flow 
including  threshold  of  sediment  motion  under  wave  action  and  combined  current 
and  wave  motion  is  given  by  Hales  (1980);  these  constitute  PARTS  IV  and  VI  of 
Report  2.  Concepts  and  ideas  that  are  required  for  understanding  siltation 
in  ocean  areas  are  discussed  below. 

15.  Boundary  shear  stress  in  sediment  transport  problems  is  of  primary 
importance.  Hales  (1980)  states  that  oscillatory  fluid  motion  associated  with 
surface  gravity  waves  exerts  shear  stresses  on  the  bottom  that  are  often 
several  times  larger  than  shear  stresses  produced  by  unidirectional  currents 
of  the  same  magnitude.  He  further  states  that  shear  stresses  produced  by  wave 
motion  may  put  sediments  into  suspension  where  they  can  be  transported  by 
currents  of  a  magnitude  insufficient  to  initiate  sediment  motion. 

16.  In  riverine  environments  the  threshold  for  the  initiation  of  sedi¬ 
ment  motion  has  been  reasonably  well  documented.  However,  this  Is  not  the 
case  for  sediment  movement  under  oscillatory  wind  waves.  This  is  primarily 
related  to  the  difficulty  of  taking  measurements  under  these  conditions.  Most 
researchers  in  this  area  apply  the  curves  of  Shields  for  unidirectional  flow 
to  the  wave  condition  considering  that  there  might  be  some  error  since  these 
accelerating  currents  exert  a  larger  shear  stress  than  does  a  steady  flow  of 
the  same  magnitude.  The  evidence  is  overwhelming  that  the  Shields  function 
with  its  limitations  serves  as  a  reliable  and  quite  general  criterion  for  the 
threshold  of  sediment  movement  under  water  waves. 

17.  Komar  and  Miller  (1974)  developed  relations  between  threshold  of 

grain  movement  defined  by  a  grain  diameter  d  and  density  and  wave  period 

T  and  the  orbital  velocity  at  the  near-bottom  u  .  All  the  data  used  in 
this  development  came  from  laboratory  experiments.  Another  consideration 
stressed  in  the  review  by  Hales  Is  that  the  interaction  of  wave  trains  of 
differing  periods  under  natural  conditions  may  generate  instantaneously  higher 
velocities,  and  sediment  motion  may  occur  at  lower  velocities  than  those 
implied  in  this  analysis. 

18.  Oscillatory  flow  produces  complex  flow  patterns  that  vary  continu¬ 
ally  in  magnitude  and  direction.  Linear  wave  theory  describes  this  condition 
as  one  with  no  net  transport.  However,  various  nonlinear  effects,  such  as 
wave  asymmetry,  and  wave-induced  net  transport,  modify  this  equilibrium  and 


to  easterly  waves  from  offshore  and  to  waves  generated  within  the  lower  bay 
which  may  reach  heights  of  over  4  ft,  especially  with  northerly  winds. 
Sediment  distribution 

13.  The  dominant  surficial  sediment  of  Chesapeake  Bay  entrance  Is  a 
homogeneous  gray,  fine  to  very  fine  quartzose  sand,  usually  well-sorted  and 
often  silty  (Melsburger  1972).  figure  A4  shows  the  surface  sediment  distribu 
tlon  In  the  study  area. 


Figure  A4.  Surface  sediment  distribution  In  the  study  is  based 
on  gross  characteristics  (Melsburger  1972) 


TIDAL  RANGE,  FT 


Figure  A3.  Tidal  range,  ft,  20-year  tidal  range  statistics, 

Virginia  Beach,  Virginia 

10.  Bosserman  and  Dolan  (1968)  studied  857  storms  for  the  period 
1942-1967  and  constructed  a  storm  probability  curve  for  three  hindcast  deep¬ 
water  wave-height  intervals.  During  the  Ash  Wednesday  storm  of  March  1962, 
the  maximum  computed  significant  wave  height  was  about  10  m.  The  severest 
storms  occur  from  January  through  March. 

11.  Hurricanes  generally  move  from  southwest  to  northeast  in  the  study 
area  (Ho  and  Tracey  1975).  The  number  of  hurricanes  affecting  the  area 
increases  from  Cape  Henry  to  Cape  Hatteras. 

12.  Waves  on  the  open  coast  south  of  Cape  Henry  as  measured  by  a  Coastal 
Engineering  Research  Center  (CERC)  wave  gage  at  Virginia  Beach  are  less  than 

3  ft  high  more  than  90  percent  of  the  time.  Most  of  the  bay  entrance  is  open 


Figure  A2.  Gross  morphology  of  the  bottom  in  Chesapeake  Bay  entrance 
area.  Soundings  are  in  feet  (Meisburger  1972) 

relative  to  land  at  Norfolk,  Virginia,  was  +4.4  mm/yr.  From  1940  to  1978,  the 
average  was  +3.7  mm/yr,  or  about  15  percent  less  than  the  1928-1978  average. 
Meteorological  and  wave  influences 

9.  Mean  annual  wind  velocities  are  16  km/hr  at  Cape  Henry.  Winds  at 
or  above  this  speed  are  predominantly  onshore  from  the  northeast  and  occur 
most  frequently  during  the  winter  months. 


A3 


Fluid  Dynamics  of  the  Chesapeake  Ray  Entrance 


Morphology 

3.  A  major  portion  of  the  Chesapeake  Bay  entrance  is  less  than  35  ft 
in  depth,  but  greater  depths  occur  in  the  channels  where  there  are  a  few 
isolated  closed  depth  depression  contours. 

4.  The  entrance  is  10  miles  across  and  extends  from  Fisherman  Island 
near  Cape  Charles  to  Cape  Henry;  however,  the  main  natural  inlet  channel  is 
less  than  2  n.m.  in  width.  This  channel  is  deepest  off  Cape  Henry  and  extends 
southeastward  for  about  5  n.m.  The  Tail  of  the  Horseshoe  is  a  sandy  shoal 
between  the  Thimble  Shoal  and  Chesapeake  Channels.  The  Cape  Charles  terrace 

is  characterized  by  numerous  secondary  morphological  features,  among  which 

★ 

linear  shoals  and  semiclosed  depressions  are  most  common  (Meisburger  1972  ). 
Figure  A2  shows  the  gross  morphology  of  the  bottom  in  the  Chesapeake  Bay 
entrance . 

Entrance  and  shelf  circulation 

5.  The  circulation  in  the  Chesapeake  Bay  is  primarily  a  result  of 
tidal  action  and  wind. 

6.  Net  nontidal  circulation  in  the  Chesapeake  Bay  Bight  (Cape  Henlopen, 
Delaware  to  Cape  Hatteras)  was  documented  by  Harrison,  Brehmer,  and  Stone 
(1964).  Their  bottom  drifters  released  on  the  shelf  up  to  distances  as  much 
as  approximately  40  n.m.  offshore  tend  to  drift  shoreward  with  some  even 
having  a  tendency  to  travel  toward  and  enter  the  Chesapeake  Bay. 

7.  The  tides  in  the  Chesapeake  Bay  entrance  are  semidiurnal  with  a 
mean  range  of  approximately  3  ft  and  a  spring  range  of  3.5  ft.  On  the  outer 
coast  of  Virginia,  adjacent  to  the  bay  entrance,  mean  and  spring  ranges  are 
about  3  and  4  ft,  respectively.  Although  the  magnitude  of  tidal  currents 
will  vary,  they  are  generally  in  the  range  of  1  to  2  knots  maximum.  Fig¬ 
ure  A3  is  the  distribution  of  predicted  tidal  ranges  for  Virginia  Beach 
summarized  for  the  20-year  period  from  1  January  1956  to  31  December  1975. 

8.  Sea-level  change  data  are  not  available  in  the  study  area.  However, 
tide  gage  records  from  Norfolk,  Virginia,  for  the  period  1928-1978  provide 
some  information  on  the  magnitude  of  changes  that  have  occurred  at  the  nearest 
stations  (Hicks  1981).  From  1928  to  1978,  the  average  rate  of  sea-level  rise 


gridded  by  automated  techniques.  A  section  map  of  the  area  was  also  made  that 
determined  regions  where  the  surveys  would  be  compared  (Figure  A10). 

36.  The  algorithm  calculates  a  surface  using  an  interpolative  procedure 
based  on  input  data  of  nearest  neighbors  to  the  calculation  location,  and  the 
calculation  procedure  provides  an  exact  fit  to  the  input  data  at  all  input 
data  locations. 

37.  Tables  A3  and  A4  summarize  the  calculated  information  from  each 
section.  Figure  All  shows  the  depth  change  (ft)  in  each  section  and  Figure  A12 
shows  the  contour  plots  of  the  two  surveys.  Figure  A12  does  indicate  a  trend 
of  natural  deepening  of  the  area  where  the  proposed  Atlantic  Ocean  Channel  is 
to  be  built. 

Analytic  methods 

38.  The  Moriches  Inlet  method  was  the  only  purely  analytic  method  that 
could  be  applied  to  this  study.  In  this  method,  a  transport  ratio  is  calcu¬ 
lated.  The  expression  used  compares  the  sediment  transport  potential  in  the 
dredged  cut  with  the  sediment  transport  on  the  bar  before  dredging.  For  this 
application  a  transport  ratio  of  0.8785  was  obtained  which  means  the  pro¬ 
posed  channel  would  cause  about  a  12  percent  reduction  in  sediment  transport 
capability  along  the  channel.  In  the  Lean  method,  some  terms  in  the  equation 
could  not  be  characterized.  The  Oregon  Inlet  method  could  not  be  applied 
since  some  of  the  key  input  parameters  were  not  available.  Generally,  the 
analytic  techniques  considered  require  more  information  than  is  available  for 
the  Atlantic  Ocean  Channel. 


Empirical  methods 

39.  The  simplified  shoaling  rate  method  and  the  Volume  of  Cut  approach 

were  not  applied,  since  these  techniques  require  an  existing  dredged  channel 

so  that  an  extrapolation  can  be  made.  Application  of  the  Gole  method  produced 

a  shoaling  rate  of  1.7  cm/yr  distributed  over  the  proposed  Atlantic  Ocean 
3 

Channel  or  88,000  m  /yr.  Details  of  the  calculations  for  the  Gole  method  are 

given  in  the  addendum  to  this  appendix.  Based  on  the  1980  survey  depths 

3 

(7.5  miles  of  channel),  there  is  an  average  of  750,000  m  /mile  of  dredging  re¬ 
quired  to  initially  develop  this  channel. 

A0.  The  method  developed  by  Boicourt  (1981),  which  indicated  a  threshold 
velocity  of  30  cm/sec  for  sediment  transport,  was  used  as  a  guide  to  determine 
when  sediment  would  move.  The  technique  proposed  by  Kadib  (1976)  was  not  ap¬ 
plied,  since  it  includes  parameters  that  could  not  be  quantified  for  this 


Figure  All.  Section  data  (depth  change  between  surveys,  ft) 


study.  Applying  the  general  concepts  proposed  by  Ludwick  (1981),  a  deposition 

3 

value  of  3.1  cm/yr,  or  163,000  m  /yr  was  determined.  Details  of  the  calcula¬ 
tions  used  in  the  Ludwick  approach  are  given  in  the  addendum. 

Summary 

41.  The  analysis  of  the  physical  model  data  indicated  that  velocities 
sufficient  to  move  fine  sand  exist  only  during  spring  tides  and  last  for  only 
2  hr  during  a  tidal  cycle.  From  the  analysis  of  the  predicted  tidal  ranges 
for  20  years  at  Virginia  Beach,  the  spring  tide  sufficient  to  move  sediment 
exists  only  approximately  10  percent  of  the  time.  Therefore  tidal  currents 
acting  alone  would  not  cause  any  significant  change  in  the  bed. 

42.  The  analysis  of  20  years  of  wave  data  indicates  that  at  the  present 
average  depth  of  the  channel  (56.96  ft)  wave  orbital  velocities  are  greater 
than  30  cm/sec  9.3  percent  of  the  time.  For  the  60-ft  Atlantic  Ocean  Channel, 
the  wave  orbital  velocities  are  greater  than  30  cm/sec  8.39  percent  of  the 
time.  Therefore  the  ability  of  the  waves  to  place  fine  sand  into  suspension 
exists  only  during  a  small  portion  of  the  year,  and  the  deepened  channel 
would  reduce  that  duration. 

43.  One  of  the  most  important  analyses  for  this  study  was  the  hydro- 
graphic  survey  comparison.  In  the  area  of  the  proposed  channel  there  was 

a  long-term  scour.  The  real  significance  of  this  comparison  is  to  demonstrate 
the  relative  stability  of  the  area. 

Conclusions 

44.  As  a  result  of  the  study,  it  is  concluded  that: 

a_.  Based  on  the  hydrographic  survey  analysis,  the  site  of  the 
Atlantic  Ocean  Channel  is  stable  to  slightly  deepening  at  a 
rate  of  about  3  cm/yr. 

_b.  The  project  depth  along  the  proposed  channel  alignment  is  at 
most  10  percent  greater  than  the  natural  depth.  The  increased 
depth  will  initially  cause  deposition;  and  the  channel,  once 
dredged,  will  have  an  equilibrium  depth  somewhere  between  the 
dredged  depth  and  its  present  depth. 

£.  The  wave  climate  is  mild.  The  majority  of  the  waves  do  not  af¬ 
fect  bottom  sediments.  This  condition  will  not  be  appreciably 
affected  by  the  deepening. 


cL  The  average  rate  of  infill  should  be  small,  ranging  somewhere 
between  2  cm/yr  (Gole)  to  3  cm/yr  (Ludwick).  These  rates  of 
infill  correspond  to  an  annual  maintenance  dredging  requirement 
in  the  range  of  about  115,000  to  200,000  cu  yd.  For  planning 
purposes,  the  higher  rate  is  recommended. 


Chesapeake  Bay  Entrance  20-Year  Wave  Summary  for  All  Directions,  Station  10 


Shoreline  angle  =  28.0  deg  azimuth.  Water  depth  =  10.00  m.  Percent  occurrence 
or  all  directions,  hs  =  significant  wave  height. 


Virginia  Reach  20-Year  Wave  Summary,  for  All  Directions,  Station  77 


reline  angle  =  342.0  deg  azimuth.  Water  depth  =  10.00  m.  Percent  occurrence  (X100)  of  height  and  period 
all  directions,  hs  =  significant  wave  height. 


Section  Data  Summar 


Section  Data  Summary 
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Volume  Change,  yd  ,  and  Depth  Change  Between  Surveys,  ft 


ADDENDUM:  ATLANTIC  OCEAN  CHANNEL  CALCULATIONS 


Gole  Method 


1.  The  Gole  method  is  an  empirical  method  developed  for  the  prediction 
of  siltation  in  harbor  basins  and  approach  channels.  The  method  is  designed 
to  treat  suspended  loads  and  is  based  on  the  analysis  of  prototype  and  model 
data  and  theoretical  studies. 

2.  Included  in  the  assumptions  made  for  developing  this  method  are  the 
following : 

a.  The  suspended  sediment  capacity  of  current  is  proportional  to 
the  square  of  velocity. 

b.  The  mechanism  of  turbulence  that  really  keeps  silt  in  suspension 
is  not  considered. 

c.  Flow  is  assumed  to  be  perpendicular  to  channel  alignment. 

(Based  on  the  physical  model  observations,  the  flow  in  the 
vicinity  of  the  proposed  channel  tends  to  be  parallel  to  the 
channel.  However,  there  are  sufficient  periods  of  time  during 
the  tidal  cycle  in  which  flow  across  the  channel  exists  to 
justify  the  application  of  this  technique.) 

3.  The  equation  is  as  follows: 


S  = 


K  L  v  C  t  d 


where 


S  =  effective  silt  load  ^ML/T^J 

K  =  empirical  coefficient  (use  0.31  for  approach  channels) 

L  =  channel  length  M 

v  =  water  velocity  approaching  channel  ^L/T^ 

C  =  silt  concentration  in  water  column  [m/l3]  approaching  channel 
t  =  time  [T] 

d  =  depth  of  flow  approaching  channel 
B  =  width  of  the  channel 
Vq  =  particle  fall  velocity  [L/T] 

V  =  water  velocity  across  channel 
D  =  channel  depth  H 


A2A 


4.  The  specific  values  used  to  make  the  estimate  for  siltation  were: 

L  =  17,059  m 

v  =  0.46  m/sec  (average  representative  value) 

C  =  0.007  kg/m3 

t  =  3.16  x  107  sec  (1  year) 

d  =  15.90  m  (52.15  ft) 

V  =  0.40  m/sec 

D  =  18.29  m  (60  ft) 

B  =  304.8  m  (1,000  ft) 

V  =  0.001  m/sec 
o 

Using  these  values 

S  =  100  x  10^  kg  annually 

Assuming  that  the  density  of  the  deposited  material  is  1134  kg/m3,  the  volume 

3 

of  siltation  is  88,000  m  annually,  or  an  average  over  the  channel  of 
1.7  cm/yr. 

Ludwick  Method 


5.  The  Ludwick  method  is  based  on  the  expression  that  approximates  the 
deposition  of  fine-grained  sediment  by  sett  Ling  in  the  presence  of  a  moving 
current 


dm  -  -  /. 

-  =  c  w  !  - 


Lo 


Lc 


For  deposition 


wliere 

m  =  mass  of  sediment  deposited  per  unit  area  of  bed  [m/L2] 
t  =  time  [l] 

c  =  deptL)-mean  concentration  of  suspended  sediment  [l/t] 
w  =  mean  settling  velocity  of  suspended  sediment  [l/t] 

Lo  =  shear  stress  of  bed,  due  to  moving  sediment  and  water  [f/l2] 

Lc  =  critical  shear  stress  above  which  sediment  capnot  be  deposited  [V/L"J 

6.  If  it  is  assumed  that  wlien  Lo  <  Lc  ,  Lo/Lc  =  0  ,  then  integrating 
yields  m  =  c  w  t  .  If  this  mass  of  deposited  sediment  is  transformed  into  a 


thickness  accumulation  rate  R  ,  where  R  =  m/pst 
then 


T  =  (c  w  /p  )  (Ludwick  1981) 
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c  =  depth-mean  concentration  of  suspended  sediment  in  the  water  [M/l  ] 

w  =  mean  settling  velocity  of  suspended  sediment  [ L/T  ] 
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Pg  =  in-place  bulk  density  of  sediment  [M/L  ] 

Values  of  constants  (Ludwick  1981): 
c  =  7  x  10  g/cm 
w  =  0 . 1  cm/sec 
p  =  1.134  g/cm^ 

R  =  7  x  10  g/cm  x  0.1  cm/sec/1.134  g/cni 
R  =  6.17  E-07  cm/sec 
R  =  2.22  E-03  cm/hr 

Assume  deposition  occurs  only  during  the  time  of  slack  water  for  only  1  hr. 


365-1/4  days  x  24  hr/da^ 
6.21  hr 


=  1,411.59  slack  waters  per  year 


2.22  E-03  cm/hr  x  1,411.59  hr  =  3.1  cm/yr 
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This  results  in  a  prediction  of  163,000  m  /yr. 


